Integral Representation of Zero-Memory 
Nonlinear Functions 
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(Manuscript received July 30, 1962) 

Integral representation of zero-memory nonlinear functions offers prom- 
ise as an analytical method for nonlinear control systems study. A review of 
work performed at Bell Laboratories and elsewhere on the use of these 
representations is presented, with particular emphasis on nonlinearities 
often encountered in feedback control systems. In general, the integral 
representations are useful only insofar as the resulting expression can be 
readily evaluated. The use of Bennett functions systematized the formulation 
of these integrals. The numerical results of a large class of the integrals 
can then be given by the tabulated Bennett functions. A comprehensive 
bibliography is appended. 

I. INTRODUCTION 

Integral representation of zero-memory nonlinear functions has been 
extensively used by Bennett, Rice and others (see References) in the 
solving of problems such as the finding of modulation products when 
one or more sinusoids appear at the input, and the finding of the output 
autocovariancc function when sine wave and random noise are applied. 
In relation to the necessary calculations which occur in the use of these 
integral representations, a class of functions known as Bennett functions, 
after W. It. Bennett, has been defined. A selected representation of these 
functions has been tabulated and plotted. 

While the original studies were carried out in relation to problems 
encountered in communications, the methods and the results can cer- 
tainly be applied to advantage in control problems. Some work in this 
regard has been done by J. C. Lozier in unpublished notes on the analysis 
of the oscillating control servomeehanism. On the whole, however, it 
appears that these approaches are not known to investigators in the 
controls field. The present paper represents an attempt to summarize 
in a unified manner the work that has been done and to indicate the 
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scope of applications and limitations of the integral representations, 
particularly with respect to controls usage. 

II. INTEGRAL REPRESENTATION ARISING FROM FOURIER TRANSFORMS 

It is known that the function 

r / \ u . u r sin uX , i « \ 

Z TT J A 

is discontinuous in its first derivative with respect to u, its value as a 
function of u being: 

/i(u) = u u > 

= w < 0. 

The plot of /i(w) vs u is in the form of an ideal half -wave rectifier. 

Using fi(u) as a basic unit, other discontinuous functions can be 
generated. For example 

. / \ 1 . 1 /*°° sin wX ,. * . n 

Mu) - - + - / — r— d\ = 1 u > 

Z 7T •'o A 

= u < 

and is in the form of an off-on relay as a function of u. 

From (3) the bang-bang type of relay is readily created as: 

. . , 2 A f 00 sin uX . . n 

/»(«) = — —7- = 4 w > 

7T J A 

= -A m < 0. 

A relay with dead zone is 

. , N A /"" sin (w - c)X + sin (u + c)X , 
U(u) = - \ r rtA 

T J A 

= 2A rsmu\cosc\ dx= Q _ c<u< c (5) 
* J ° x = -A u < -c. 

A limiter (linear characteristic with saturation) is 

/,(«) = ? P sin ^ Xsin ' 4X rfX = u" -/1< I < A (6a) 
* " x " = /l 1* > A 

= 2u r cos uX sin AX ^ + 1A T sin 11X cos AX dx (fib) 

IT Jo X T Jfl X 
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Equation (6b) is readily obtained by manipulating two functions of 
the form of (1). That (6a) is equivalent to (6b) is seen by integrating 
(6a) once by parts. Other discontinuous functions can be generated 
from the above five functions by appropriate shifting (bias) of each 
individual characteristic, or by combining several characteristics. In 
fact, simply multiplying by an appropriate g{u) can create quite general 
discontinuous characteristics. 

It is noted that u may be viewed as the input to the nonlinear element, 
and f(u) then gives the response to this input. If u{t) is a function of 
time, for each u(t\) the function f[u(ti)] yields the instantaneous value 
of the output (i.e., f(u) is a functional of u). While f(u) is no more 
convenient for use in the evaluation of the output as a function of time 
than equations of the form (2), giving the discontinuous function as a 
set of equations, it is very useful for the purpose of spectral analysis 
since f[u(t)] is in a compact form suitable for Fourier series expansion. 

As an example, we seek to find the output spectral component for the 
relay with dead zone, when input is in the form u = P cos x. 

Using (5), the output Fourier coefficients are found by: 



T J—r T Ji) 



sin (XP cos x) cos cX 



d\ > cos nx dx 



2 A f°d\ 



= 2A f 

•7T 2 in 



7T 2 Jo X 




cos cX 



L sin ( 



XP cos .r) cos nx dx 



n even 



A. (-i)<-" f J^P\ °°s <* (/x „ odd 

47T •'u A 



where ./„(z) is Bessel function of the first kind of order n. 
Since it is known that 1 



1 j • -i fb 

- cos < n sin 

n \a 



f J n (al) cos /;/. „ 

/ : dt = \ n nir 

J u t a cos — 



then, for n odd: 
A 



«„ = 4 -(-d 

= 



n\b+ Vb°- - a 2 }" 
"■{i«*[nam- , (£)] 



b < a 
b > a, 

c< P 
P <c. 



(7) 
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Certain of the nonlinear characteristics expressed in the integral form 
are also amenable to a double frequency type of analysis in which the 
input is of the form 

u(t) = P cos (toi< + 0i) + Q cos («2< + 6 2 ). 

Bennett 2 in particular has contributed extensively to double frequency 
studies. 

In control systems analysis, a double frequency study becomes neces- 
sary in (a) the oscillating servomechanisms 3 and (b) the dual input 
describing function approach to closed loop servos. 46,6 In what follows, 
the fundamental components (i.e., components in coi and w 2 of the out- 
put) from a bang-bang type of relay are found. The approach follows 
closely that of Lozier. 

The input u is a function of two frequencies wi and o> 2 ; this is brought 
to light by setting x = tait + ft , y = o> 2 t + 2 , and letting Q/P = k. 

Thus 

u(x, y) = P(cos x + k cos y). (8) 

In passing through a bang-bang type relay, it is recognized that the 
amplitude P in (8) does not influence the output; thus without loss of 
generality it may be set to unity. 

The output /(m), written asf(x,y) is: 

f(x,y) = +A cos x + k cos y > 
f(x,y) = —A cos x + k cos y < 
which may be expressed as a double Fourier series 2 ' as: 

00 00 

f(x,y) = E E [A±mn cos (mx ± ny) + B ±mn sin (mx ± ny)\ (0) 

ii =-0 »«=0 

where 

ZW~ J—r J-r 

B±mn = 7T-o [ I f( X >V) Sm ( MX ± nl j) ^ *" ^ 

2 71" J-x J-t 

Aoo = 3 f P/(*tf) * *». B «m - 0. (9c) 

7T" J-7T J-T 

From symmetry of the bang-bang relay, B± m „ = for all m and n; 
moreover the integral representation of (4) can here be used, thus: 
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a A [* [ W V f* d\ sin (cos x + k cos y) X~| 

A± "» - -? L L L-'o 1 J 

•cos (mx ± ny) dy dx. 

The interchange of integration can be carried out here in view of the 
finite limits of the outer integral and the bounded nature of the inner 
integral, whence: 

A ±mn = -% I — I dx I dy {sin [(cos a; + k cosy)\] cos (mx ± ny) I. 

7T' J A J~r J-k 

Upon expanding, collecting nonzero terms, and integrating, the result 
is (as A +mn = A_ mn for all m, n, the ± signs are henceforth dropped.) : 

A mn = M (.d^-d/^ r jn(k\)j m (\) m n odd 

7r ''o A (10; 

= otherwise 

where use has been made of the following definite integrals: 8 ' 1 

2 f* 12 

- I cos (z sin tp) cos 2rupd<p 

IT Jo 

„ 2 r' 8 <Ua) 

= (—1)" - / COS (Z COS ip) COS (2n<£>) rfcp = J<> n (z) 

7T "Ml 

2 Z"' 2 

- / sin (z sin <p) sin (2n + l)<pd(p 
ir Jo 

.2 /•'" (Ub) 

= (— l) n '- / sin (z cos (p) cos (2n + l)<pd<p = J2n+i(^)- 

7T ''O 



The integral in (10) may be evaluated by means of formulas attributed 
to Sonine and Schafheitlin' (also known as the Weber-Schafheitlin in- 
tegrals 8 ), the result being expressed in the form of hypergeometric func- 
tions F(a,(3,c,x) : 

„ (n + m — r + 1 

f K Jn(a\)JJ b\) ^ _ al \ 



d\ = 



2^+. r( -» + m 2 + r+ > (». + !) (12a) 



,, (n + m — r + 1 n — m — r ■+• 1 /V 

if n + m - r + 1 >(),/•> - 1 , and < a < b, 
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r c ° «/ w (aX) t / OT (aX) dx 

gfr(r)r (» + 'V + 1 ) (12b) 



2r / - n + m + r + l \ p / n + m + r+l \ p / n - to + r + 1 \ 

if n + w + 1 > 0, r > 0, and a real; and 

FT (n±JH—I±l) 
f K Jn(a\)JM\) V 2 J 

h " xr «~»r ( n ~ m 2 +r + 1 ) r(« + i) (12c) 



( n + to - r + 1 -ft + to - r + 



1 j_ i 6 ^ 



if (n + to - r + 1) > 0, r > -1, and < b < a. 
Accordingly, 

, __ 2 A / 1 v(n+m-l)/2 \ 2 / 



_ /2 — n + to\ „, . ., N 

r ( 2 ) r < B + 1 ) (13a) 

p /n+_ro ^ r^-jn ^ b + ^ /c A fQr fc < x 



_ 2A , - •, ( n + m _D/2 

7T 



(n. + to\ 
2 ; 



(13b) 



r ( 2 " w 2 + m ) r ( 2 + n 2 +m ) r (^f^- 2 ) 



for fc = 1 



_ ^- / _ 1 \ (n+n-1) - \ * / 



2 A / .. x („ +n _D/2 



«(^±_ 2 )r( m + 1 ) (13c) 



_ (n 4- to to — ft . . / 1\ \ f ,. «. , 

,/? ' Hr" ' -2~ > m + 1 '\k)) for /c > L 
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The three cases of (13) are essentially equivalent if one recognizes that 

r(c)r(c - a - b) 



F{a,b,c,\) = 



T(c - a)T(c - b) 



and that, for k > 1, the situation is identical to that of inverting the 
role of n and m, and defining a new quantity A;' = 1/k. 

The output fundamental components are i4i and A m ; from (13) one 
has: 



Ak = m r jk Fl t 



*r(|)r(i)-6'~ 2i,M ' ! ) 



(14a) 



_ 4 A f 1 2 > , ■! ■ > 
"Tl'-i' "Si* -256* 



3,4 5 , 



A 0l = 






(14b) 



Considered together with the input (8), this yields Lozier's oft-quoted 
result," 10 that the equivalent "gain" of the relay, for small values of 
k, is 6 db higher for the "carrier" than for the "signal." 

It is not difficult to see that the Weber-Schafheitlin integrals also 
occur for two frequency inputs applied to the characteristics (1) and 
(3), but that integrals of the form 



f 

J a 



M«K)J.(bX) fsin .4X1 dx (15) 



fsin AX) 
\cos A\) 



occur for characteristics of (5) and (6). Moreover, inputs with more 
than two frequency components will result in generalized Weber-Schaf- 
heitlin integrals of the form 

\ HJi(ai\)\- r d\ (16) 



t Equations (14a) and (14b) can also be expressed in terms of the complete 
elliptic integrals for which tables are available. See Refs. 3, 21, and 22. 
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for characteristics (1), (3), and (4), and generalized integrals of the 
form 

/"n^rf n ^l& (it) 

Jo £=i [coSilXJ 

result for characteristics (5) and (6) . 

Unfortunately no general solutions have been found to represent 
(15), (16) or (17) in known functions. In such cases, numerical solutions 
can be used. Numerical solutions in terms of Bennett functions and their 
tabulation are described in Section IV of this paper. 

III. INTEGRAL REPRESENTATION ARISING FROM LAPLACE TRANSFORM 11 - 12 

The integral representation of Section II is closely related to the 
Fourier transform. An alternate approach using the Laplace transform 
is more convenient in some cases and has been extensively used by 
Bennett and Rice, among others. We mention some results of this ap- 
proach for the sake of completeness. 

Expressing the output of a nonlinear device in response to an input 
u as/(w), it is possible to find the (possibly two-sided) Laplace trans- 
form of f(u), denoted F(s), or, 

F(s) = f e- au f(u)du. (18) 

J— 00 

The inverse transform is then 

f( u ) = -J-. [ e " 8 F(s)^ (19) 

where C is some suitably chosen contour of integration. If F{s) exists, 
then (19) is an explicit expression for f(u) which may be used to ad- 
vantage. In the case of solving for modulation products, f(u), written 
explicitly in x and y, may thus be used directly to compute the double 
series coefficients. 

To compute A mn , for example, using double Fourier Series expansion 
in response to an input u = P cos x + Q cos y, one has 



A - e j^ 

47T- 



f f cos mx cos ny dx dy [-L. / e' (P coa * +Q C0B v) F(s) *]* 

= ^n_ f F ( s ) ds f g" co. x cQS mx dx f e »Q cos „ cQS ny dy 

8ir 3 j ->c J-t J -r 



t « m is the Neumann factor, denned as: 

e m = 1 m = 
e m = 2 m = l,2, 
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Since 



.' — n 



•K •'O 



[ e* COBa cos na da = J n (z) (20) 

and letting s = ju, one obtains, 

A mn = ^y- + - f F(j u )J m (P u )J H (Qu) dm (21) 

and the required coefficients are evaluated by contour integration. 

The above result can readily be generalized. For example, where 
there is a dc bias of b units superimposed on the P cos x + Q cos y in 
the input, (i.e., u = b + P cos x + Q cos ?/), the net result is to insert 
a factor e 1 w under the integral of (21). 

Inputs of the form u = b + ]0-i ^*> cos - T » w iN result in coefficients 
of the form 

r r 

II e»i i exp 2 ». r 

^n,,..n r - — ^— / C*-f(i«) n^n^P.'O,) do) (22) 

7T •'c i = l 

whenever /(it) is Laplace transformable. The contour C is a function 
only of the nonlinear device, as may be expected. The Laplace transform 
of several ordinarily encountered ideal nonlinear devices, as well as their 
associated contour of integration C has been given in Rice 11 in his ap- 
pendix 4 A. 

The nonlinear devices expressed as in (19) may be used fruitfully in 
certain investigation in noise problems. These are briefly described here. 
Reference may be made to Rice's classic papers of 1944 and 1945. n 

For inputs that include narrowband noise, the input waveform will be 
of form 

u = R cos (m m t + 0) R ^ 

where R and 6 are functions of time whose variation is slow as compared 
to cos u„,t. (lo„,/2t is approximately the midband frequency.) 
The output /(w) then is 

/(w) = ?r / W«) ex P U<» R cos (o) m t + e) dot], 

Zt j c 
By means of the relation 

c *co., _ f] e ,j" cos n^J n ( 2 ) (23) 
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the equation above may be written: 

f(u) = E MR) cos (nuj + nd) (24) 

»— o 

where 

A n (R) = «»£- ( F(J")Jn(uR) do>. (25) 

Zir J c 

In this representation, important conclusions may be reached concern- 
ing the properties of the output without undertaking laborious compu- 
tations. For A H (R) whose variation is of the order of that of R, the 
output spectrum has bands which are centered at f m , 2/ w • • • . A narrow- 
band filter centered about nf m will then yield a slowly varying cosine 
wave with envelope A n (R). A narrow-band low-pass filter will yield the 

level A (R). 

In some cases the probability density function P(R) of R is known. 
(For narrow-band Gaussian noise, for example, P(R) is the Rayleigh 
distribution.) The probability density of the output envelope A,,(R) 
is simply: 

P[A n (R)]-"i^r. (26) 



dAJR] 



dR 



Another application in which the representation of (19) is useful is 
the calculation of the autocovariance function of the output of a zero- 
memory nonlinear device. From this the output power spectrum is found 
by taking the Fourier cosine transform. 

The autocovariance function of the output is : 

¥(r) - lim I f T f[u(t))f[u(t + r)] dt . (27) 

r-»» I Jo 

By (20): 

¥(t) = lim 7-4s [ [ F(jm)exp[jmu(t))du:i / F(M 

T-oo 47T 2 7 -'o -'c J C 

•exp [jw*u(t -f- t)] dioi dt. 
If an exchange of limits is justifiable, the above becomes 

¥(t) = 7—, / F(jui) dm / F(jui) don 

47T" J C J C 

• lim = / exp [juiu(t) + jwtu(t + r)] dt . 
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The quantity in the bracket is the (time) average value of 

expj[wiu(t) + Uiu(t + t)] 

which, in the event that u(t) satisfies the ergodic hypothesis, is equal to 
the characteristic function of the two variables u(t) and u(t + t). De- 
noting this quantity by g(ui , &* , t), one has: 



*(r) = j-„ / F(M dun I FO'«2)ff(coi,a- 2 ,T) dut. 



(28) 



This gives an interesting approach to the computation of the output 
autocovariance function. 

It is interesting to note, incidentally, that the characteristic function 
of u(t) = P cos pt is 

Jo\P\/(jil 2 + OJ 2 2 + 2a)lC0 2 COS J)t) , 

and for 

u{t) = P cos pt + Q cos §£ 

where p and 5 are incommensurable, the characteristic function is 

Jo(PVcoi 2 + W2 2 + 2«iw 2 COS pr) X Jo{QVo)i* + o>2 2 + 2wico 2 COS gr). 

Here, as elsewhere, one is limited by his ability to integrate. The auto- 
covariance function, however, has been solved for particular nonlinear 
characteristics, for example, the square-law device. 

IV. NUMERICAL SOLUTIONS AND BENNETT FUNCTIONS 

Since it has not yet been found possible to express the modulation 
coefficients in a more general case in terms of known functions, it is 
often necessary to resort to numerical computations. The numerical 
approaches have been tackled by Sternberg, Kaufman, Feuerstein, 
Shipman, among others. 1319 Some of their results have been tabulated 
and a class of generalized functions encoimtered in these investigations 
are christened Bennett functions. 14 

The original approach of Sternberg and Kaufman is along the lines of 
direct integration, summarized below. 

If the output f(u) can be expressed in the form of a continuous N + 1 - 
sided polygonal function over a closed interval — a ^ u ^ a, i.e., 

f(u) = f(-a) + £ (jiU-oiu - Ui ) i = 1, ... N (29) 
1=1 

where m, and g t are constants, u t being the "break-points" of the polygo- 
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nal function 

— a ^ Mi < Ms < Ma ■ ■ • Us < a (30) 

and U-t(u — w.) are unit ramp functions: 

U-i(u - Mi) =0 M < Ui 

U-t(u - Ui) = M - ?*,- M ^ M,- (t = 1, 2 • • • 2V> 

If the input is of form w = P cos a; + Q cos //, one can confine his at- 
tention piecewise to N functions of the type 

fi(x,y) = f(P cos x + Qcoswjm,-) i = 1, 2 ■ • • N. (32) 
The over-all function is then 

N 

f(x,y) = f(-a) + £ 9/<(Ssir). (33) 



Factoring out P in each term, and introducing parameters hi = Ui/P, 
fc = Q/p, we express fi(x,y) as the double Fourier series: 

oo t 

fi(x,y) = iPAao(M) + P E ^±m»(M) cos (ma; ± nw)t (34) 

>n,n=0 

where 

A±«i(M0 = — ^ f f /»(.i-,//) cos (mx ± ray) dz dy, 

m,n = 0, 1, 2 • • • ; i - 1, 2 • • • AT 

and for f(x,y) we carry out another expansion: 

/(*,!/) = hCoo + £ C±™ cos (mi ± ny). (35) 

m,n=0 

The C's and the A's are then related by 

hCoo = /(-a) + \P £ M»(M) (36) 

i =i 

AT 

c ±m „ = p EM±-(A*.*)- 



As Ah»»(M) = A_ mn (fc,-,fc) for all m and n, the ± sign can be 
dropped. 

t Z denotes a summation without the Aoo term; in addition, terms whose in- 
dex is such that m-n = are to be weighed by a factor of \. 
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By considering the function 

JP(co8X + k cos y — hi); cos x + k cos y ^ /i, 



ftfclf) 

1,0 cos x + « cos ?/ < hi 

i = 1,2 ••• N 

the zones over which integration for the evaluation of A mn must be 
carried out is seen to be bounded by the curve 

cos X + k cos y = h; h = ft,- 

over the closed square 

— ir^x^ir 



1U ■ , 

— 7T ^ 7/ ^ 7T. 

Five cases need to be considered; two are degenerate: 
(dl) 1 + A^ ^ h 

(d2) -(A- + 1) ^ /*. 

In the first instance the integrand vanishes everywhere except possibly 
over a set of zero measure, and hence the coefficients are identically 
zero. In the second instance the integration is to be carried out through- 
out the zone (excepting possibly a set of zero measure), which means 
the output is the same as the input except for a constant multiplying 
factor. 

The three nondegenerate cases are: 

(i) // < 1 + k, or 

h > 1 - k. 

The integral here is to be carried out over a zone R of the x,y plane 
bounded by a closed curve lying wholly within R . 

(ii) // ^ k — 1, or 

ft g 1 - k 

The integral here is to be carried out over a zone R bounded by two 
open curves (i.e., two opposite segments of the boundary of R also 
constitute the boundary of 7? ). 

(in) h < k — 1 , or 

h > -(*+ 1) 

The area of integration is bounded by four open curves. The integra- 
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tion can thus be carried out for these nondegenerate cases by formulas 
of the type : 

A mn (h,k) = \ \ / (cos x + k cosy — h) cos mx dx cos ny dy 

+ — I \ I (cos x + k cos y — h) cos mx dx cos ny dy 

m,n — 0, 1, 2 • ■ • 
A mn (h,k) = — / / (cos x + k cos y — h) cos ny dy cos mx dx 

-| / / (cos x + k cos y — h) cos ny dy cos mx dx 

m,n = 0,1,2 ••• 

where Ri , R-t , R3 , R< are zones appropriate for each of the cases. 

It is seen that the inner integrals can be performed, after which suit- 
able manipulation will yield a set of recurrence relationship first derived 
by Rice. 20 Except for misprints, they are : 

(to - 7i + 3)iW.»-i = -(to + n - 3)A»_i,»-i + 2mhA m , n - 1 

— 2mkA mn m,n ^ 1 
(m + n + \)A mn = -(to - n - 3)i4 m _ 2 , n - 2(to - \)kA m ^, n ^ 

+ 2(to - l)/i4 m -,. n to = 2,w = 1 

(n + m + lMr.,7. = — (n — to — 3).4 m ,„_ 2 

A (37) 

- 2(TO - 1) r ^n,-l.n-l + 2(» - 1) £ 4 m ,„-l 

TO ^ 1,71 = 2 

(n - to + 3).4 m _i, B+ i ■ -(« + to - 3)i4„,_i,„_i 

+ 2n 7 i4 m _i >B - 2n 7 4 mn to,w ^ 1 . 

With the aid of these relationships, the higher-order coefficients can 
be expressed in terms of the first four coefficients A m (h,k), A l0 (h,k), 
A i(h,k) and A n (h,k). 

For cases such as the ideal limiter, the antisymmetric condition 
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/( — ") = —/(«■) is observed; m are now symmetric and the "gains" 
0, are antisymmetric. Here: 

A m (-h,k)=A m (h,k) + 2h 

A i0 (-h,k) ss 1 - A M (h,k) 

(38) 
A 0i (-h,k) = k - A 0l (h,k) 

A mn (-h,k) = (-\y" +n A mu (h,k) (m + n > 1). 

The function A mn (h,k) are called by Sternberg the Bennett functions 
of multiplicity two and order m,n. In part II of Sternberg's paper, 14 
the functions A m (h,k), A lQ (h,k), A 0l (h,k) and A u (h,k), have been tabu- 
lated for h between — 2 and +2 in 0.2 steps and k with values of 0.001 , 
0.01, 0.1 and 1.0. The values A 2O (M0, A m (h,k), A 30 (h,k), A 21 (h,k). 
A u (h,k), A 03 (h,k) are tabulated for k of 0.1 and 1.0. All values are tabu- 
lated to six decimal places. The accuracy of the first set of tables is held 
to be to one unit in the last place, while for the second set the accuracy 
is about three units in the last place. 

The above approach is extendable to devices with continuous and 
smooth characteristics if it can be approximated in a piece-wise linear 
form. As long as the characteristic may be approximated to within a 
pre-chosen e > uniformly on the interval —a^u^a by 

Nit) 

S(u,t) = f(-a) + £ QiU-*{u - «,■), (39) 

i=i 

Sternberg and Kaufman show that the approximate modulation product 
amplitudes computed as per (33), (34) and (35) will not differ from 
the true values by more than 4e/7r in all cases, and the output will be 
within t of the true value for all time if it is obtained by summing over 
the approximate expansions. 
For a symmetrical ideal limiter 
( 

f(u) = \ gu - » < u ^ ?/„ < » < 2P,g > 0. 



guo 


U ^ -u n 


gu 


— llo < U S Uu 


gito 


U ^ Wo 



The approaches described above can be applied to the range 

— a^x^a 

a ^ 2P. 

Sternberg, in part II of his paper, 14 gives the results relating the co- 
efficients C mn and A mn as: 
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C ±mn = m + n = 0, 2, 4 • • • (40a) 

Cio - /'(/[I - 2A M (M)] ( 40b ) 

Co, = Pg[k - 2A m (ft,fc)] ( 40c ) 

C ±m „ = -2PgA mn (h,k) m + n = 3, 5, 7, • • •■ (40d) 

Here tf = Mo/P, fc = Q/P and input are in the form 

u(t) = P cos (pt + Op) + Q cos ($ + 0?) 

0<P<P+Q£2P. 

In Ref. 17, Bennett functions of the vth kind, denoted A mn ' (fc), 
are defined. These are the coefficients for the output of a vth law recti- 
fying function in response to a two-frequency input, v is usually taken 
to be an integer. A mB w (fc) for v = 1,2 have been tabulated. 17 Bennett 
functions of a given kind can be obtained from those of the lower kinds 
by means of recursion formulas. 

By extending the above, Bennett functions with multiplicities of 
three or higher (i.e., modulation coefficients when the input has three 
or more distinct frequency components) can readily be defined. For 
input of the form 

u(t) = P(cos x + h cos y + k 2 cos z) (41) 

for example, the output from a piece-wise linear nonlinear element can 
be expressed in terms of the Bennett functions of the first kind, 

A mn i(h, h ,fo), h = -p 

where, as before, u is the breakpoint for an individual segment. Simi- 
larly, for a i»th law rectifier subjected to inputs of the form (41), Ben- 
nett functions of the vth kind 

A mnl { '\h,h,h) 

can be defined. 

A number of interesting relations have been derived for the three 
frequency Bennett functions of the vth kind.' 9 These include recurrence 
relations and integrals Unking three-frequency Bennett functions with 
two-frequency ones. No tabulation of the three or more frequency 
Bennett functions is known to have been attempted. 

Relationships between the "Fourier" representation and the "La- 
place" representations for nonlinear characteristics have also been 
revealed by Feuerstein. 18 He has shown that, in many cases, the contour 
integration in the "Laplace" formulation can be reduced to integrals 
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over the infinite real line. The results are generalized Weber-Schafheit- 

lin integrals of the form ( 10) and (17). This is perhaps not a surprising 
result from intuitive grounds. 

It is interesting to note, however, that for a v-law rectifier, the Ben- 
nett function of the vth kind of arbitrary multiplicity is given by 

Aj v \... mN (h,k lf • •-,*») - v\f~'~' - f y- (v+1) cos\h 

•II./ m .(/wX)r/X 

»"=0 

for v integer, (42a) 

M ^ v + 1, and 
M + v odd, 
and 

A mo ( '\... mN (h,k u ...,*,) = v\j M —* 2 - f X- ( ' +1) 

IT J 

N 

• sin \h n «A» , ( fc<X ) d\ , ^ . 

for j* integer, 
il/ ^ v, and 
M + v even, 

N 

where M = zl m t an d fco = 1 ^ /c,- . 

t'=0 

By these formulas, the generalized integrals of (17) are related 
directly to Bennett functions. 

Feuerstein in fact did not stop with the considerations of integer v. 
The formulas of the Bennett functions of the »»th kind, with noninteger 
v and with M taking in values other than those shown above, are re- 
lated to the generalized integrals of (17), though in a more compli- 
cated form. 
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